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Abstract 

Relations between Hamiltonian mechanics and quantum mechanics are studied. It is stressed 
that classical mechanics possesses all the specific features of quantum theory: operators, complex 
variables, probabilities (in case of ergodic systems). The Planck constant and the Fock space 
appear after putting a dynamical system in a thermal bath. For harmonic oscillator in a thermal 
bath, the probability amplitudes can be identified with the complex valued phase functions f(q+ip) 
describing small deviations from the equilibrium state, when the time of relaxation is large. A 
chain of such oscillators models both the one- dimensional space (or string), and one-dimensional 
(Sj \ quantum field theory. 

en ; 

o^ ; 1 Introduction 

At the present time it becomes clear that to uncover the nature of quantum description one has 
Q \ to turn to the Planck scales Zp ~ 1.6 • 10~ 33 cm [1,2]. For example, 't Hooft came to conclusion 

that to construct self-consistent quantum gravity one should 1) admit discrete space-time and 
physical variables, 2) postulate at the Planck scales deterministic classical mechanics (CM), 3) 
admit (locally) dissipative processes [1]. Of course, at the atomic level one cannot obtain quantum 
mechanics (QM) from the classical one. It may be possible only in the models where matter, 
mechanics and space appear simultaneously [2]. For the recent activity in this direction see [3,4]. 
The problem of the Planck constant (what is the nature of hi) was not discussed in [1,3,4]. 

The aim of the article is to show that evolution of non-equilibrium states of a harmonic 
oscillator in a thermal bath is described by probability amplitudes (if the deviation from the 
S^ equilibrium is small, and the relaxation time is large). It is the heat bath that is responsible 

for appearance of probabilities and the Planck constant h. A chain of such oscillators models 
one-dimensional quantum field theory or quantized bosonic string. One-particle excitations of the 
field play the role of particles and are described by wave functions. The distances between the 
oscillators are supposed to be of order /p. Thus, ID relativistic quantum world can be built in the 
framework of non-relativistic classical mechanics. The harmonic oscillator plays the prominent 
role in this model (as is well known fields are ordered sets of harmonic oscillators). 

In section 2 the attention is drawn to the fact that Hamiltonian mechanics possesses all the spe- 
cific properties of quantum theory: operators, complex valued functions, probabilities (for ergodic 
systems). Moreover, in case of complex systems there appears also a constant with dimension of 
action and the Hilbert space. 

In section 3 equilibrium and non-equilibrium distributions for a system in a thermal bath are 
discussed. Variations of canonical variables are divided into two classes. Those, preserving the 
Gibbs distribution satisfy equations analogous to the Hamiltonian equations of motion. All the 
other variations give rise to non-equilibrium states. 
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Study of harmonic oscillator in a thermal bath (section 4) shows that phase functions f(z), 
z = (q + ip)/\/2, describing non-equilibrium states, in case of large relaxation time t r compose 
the Fock space and may be identified with the probability amplitudes. 

In section 5 a chain of harmonic oscillators in a heat bath is considered. It is shown that in 
the continuum limit one obtains ID quantum field theory. It elucidates the nature of quantum 
fields, matter (excitations of the fields) and physical space (the chain of oscillators). 

In Conclusion (section 6) the main results of the paper are summarized, and the consequence 
of the proposed model for the Universe is pointed out (disappearance of coherent excitations of 
the vacuum, i.e. "disappearance of matter"). 

2 Hamiltonian mechanics and mathematical apparatus of 
quantum mechanics 

To define Hamiltonian mechanics, one has to define 1) phase space (PS) --an even-dimensional 
manifold, 2) non-degenerate closed symplectic form on PS, 3) the Hamilton function H(q,p). The 
symplectic form 



u? 



(q,p) = E^'fc^* Ad Pk (K 

fc=i 



defines the Poisson bracket for phase functions f{q,p),g(q,p) 

n ( df dg df dg\ n d(f, g) 

^ k ' \dq k dp k dp k dq k J ^ k ' d(q k ,p k )' 

Because we are going to use non-canonical transformations, the following symbol is proved to be 
useful 

M™> = & (3) 

The equation of motion reads 

f = {f,H}. (4) 

Hamiltonian mechanics contains all the essential constituents of quantum mechanics. For the 
readers convenience we give some details used in what follows. 
Operators. Equation (4) can be rewritten in the form 

f = {H, /p,q}f = H d f, (5) 

where {H, /p,q} is an operator: {H, /p,q}f = {f,H/q,p}. For simplicity, we took n = 1 and 
^(QyP) — 1- H d is an operator in CM [5]. 

Complex functions. Let us take as an example harmonic oscillator with the Hamiltonian 

1 ,p 2 „2\ W/ 2 2\ FT 1 P 2 1-1 , /^\ 

H = -{ y iq ) = -{p +q), u = J — , p = —=, q =q^m. (6) 
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The equations of motion are 

±(q,p) = u.Jx(q,p), J=( _ 1 j, J 2 = -1 (x.(q,p): X! = q,x 2 = p). (7) 

In normal coordinates z = (q + ip)/^/2, z = (q — ip)/\/2 we have 

x(*,z) =iS'x(z,z), x(f,z) = Ux(q,p), S = ( Q _ 1 J , U = —= ( 1 /J, £> + £> = 1, 
i.e. Eqs. (7) decouple, and one may use e.g. only the second equation 



the first one can be obtained by complex conjugation. Complex canonical variables z, z are useful 
in general case [6,7]. 

The imaginary unit. Notice, that 1) matrix J in (7) is transformed into the imaginary unit 
i in (8): UJU + = iS, 2) transformation q,p — > z,z is not canonical one: {z,z/q,p} = i (though 
U is a unitary matrix; for a canonical transformation q,p — > Q, P by definition {Q, P/q,p} = 1). 

Probabilities. Complex ergodic systems are stochastic ones; the time average of phase func- 
tions is equal to their average over an ensemble. Thus, in deterministic classical theories there 
appear probabilities and (for subsystems, see e.g. [8]) the Gibbs distributions. This is one of the 
most important points of our consideration. 

Hilbert space. The Gibbs distribution allows to define the Hilbert space; for complex valued 
phase functions f{q,p),g{q,p) one can introduce the inner product [5,9] 

(/, g) = Z~ l [ fge-? H ^) d n qd n p (9) 



(here /3 = 1/kT, k - - the Boltzmann constant, T - the temperature, Z - - the normalizing 
constant). It was stressed however [9], that in QM wave functions are defined in configuration (or 
momentum) space, not in PS. The connection between the Hilbert space with the inner product 
(9) and the Fock space was not noticed. The point is: arbitrary phase function f(q,p) cannot be 
a wave function because q and p do not commute. But the phase functions f(q,p) = f(z) can be 
considered (and are considered) as wave functions in the Fock space -- now canonical variables 
are z, z (see section 4). 

The Planck constant h. For finite motion the integral 



exists. Then, the constant 



-^(i^d n qd n p (10) 



z 1/n = h (ir 



has the dimension of action, and in certain cases it can be identified with the Planck constant 
h (see sections 4,5). Of course, in general case this h depends on H, /3 and is not a universal 
constant. But in a world made of identical elements, e.g. in ID quantum field theory (linear chain 
of oscillators) this constant is universal one. 



3 Equilibrium and non-equilibrium distributions 

The Gibbs distribution and the Hamiltonian equations of motion. As was already 
mentioned, classical mechanics of complex ergodic systems is stochastic and bears the Gibbs 
distribution. Problem: find variations of canonical variables preserving the Gibbs distribution. 
From 5exp(—(3H) = one gets 

—5q t + — S Pi = £ ViHSxi = VH5x = (12) 

OQi OPi ) i 

(xj = x(gj,pj), x = x.(q 1 ,p 1 , ..., q n ,p n ))- Solution of this equation is 

5x = ClVHSt, (13) 

where Cl is some In x In antisymmetric matrix, such that 

VHVLVH = 0. (14) 

Among solutions (13) there are solutions of the Hamiltonian equations of motion. In the simplest 

n 

case (fi =® J, cUk\Q,P) — 1) equations (13) are nothing but the familiar Hamiltonian equations 

. 3H . dH 

Qi = ^— , Pi = — q— ■ ( 15 ) 

opt dqi 

It follows from (12) that an arbitrary variation <5x is the sum 

5x = 5n.±_ + <5x||, (16) 

where 5x.± are given by (13). Variations <5x^ preserve the Gibbs distribution while variations 5xy 
give rise to non-equilibrium states. 

Non-equilibrium distributions. The Gibbs distribution for harmonic oscillator introduces 
measure in PS 

d(*{z,z ) = — — e~ zz / h , h = — = — 17 

in 27r poo 

(see (6), (11)). Any other measure ji p ^ /j, describes some non-equilibrium state 

dfip = (d/ip/ d/j,)d/j, = p(z, z)dfi, p(z, z) > 0. (18) 

Now, consider variation of z, z with small non-zero 5xy in (16); we have 

uzz -> uzz - P~\cz + cz) + ..., e~^~ zz -> |e C2+ -f e-^/ ft , 

and 

d^,«) -» dn f (z,z) = \f(z)\ 2 dfi(z,z), f(z) = f c (z) = e cz ; (19) 

the higher degrees of z, z in the exponential are omitted. Then, instead of p(z, z) in (18) we have 
\f(z)\ 2 , so to study evolution of a non-equilibrium state it is enough to consider evolution of f(z). 
Evolution of a single analytical function f(z) describes evolution of a non-equilibrium state; entire 
functions f(z) are elements of the Fock space [10,11]. This is the key to the problem. 



Remark. Introduction of h by Eqs. (11), (17) may look unconvincing at the first sight. Notice, 
however, that (17) can be obtained by projection of a sphere (tp, 6) on complex plane z 

\A 2 = a ln ofl D 2 • 2 /i/o ' ar S^ = ^; R 2 sm6dipAde = d^(z,z), 
p 2pR 2 sin 0/2 

where R is radius of the sphere [2]. If PS is a sphere then its area can be identified with h = 47rR 2 
- the universal constant with dimension of action for Hamiltonian dynamics on the sphere. A 
chain of such oscillators can model a bosonic string, while a network of strings models e.g. 3D 
space [2]. Hamiltonian mechanics in this space possesses the universal constant h. 

4 Harmonic oscillator in a thermal bath and quantum me- 
chanics 

It is easy to show now that for harmonic oscillator evolution of small (in the sense of section 3) 
deviations from the equilibrium state is described by quantum mechanics if time of relaxation is 
large (t r > to^ 1 ). 

Evolution of non-equilibrium states. Hamiltonian (6) in complex variables takes the form 

'-i(;)(i!)(;)-?('M*)-?(*)(!i)(*)-> + - ) - 

(20) 
and the equation of motion for f(z) reads 

/ = {f,H/q,p} = t{f,H/z,z} = -i^f(z) = -iuJL-Hz) = H d f (21) 

or, in terms of <p = arg z 

(d t + ud lfi )f(<p,t) = (22) 

(d\z\/dt = 0). Functions f(ip—ut) solve equation (22); they describe waves on the circle \z\ = const 
and should be periodic functions. So, f(z) can be developed into the Laurent series. But in general 
case such functions are singular at z = and do not suit for description of any sensible state of 
harmonic oscillator. The proper functions are those given by the Maclaurin series 



f(z) = J2 CnZ n (z), Z n (z) = z n /Vn~\. (23) 

n=0 

Functions Z n satisfy equations of motion 

Z n = —iujnZ n , (24) 

and evidently correspond to motion with frequencies um. 

Probability amplitudes. Solutions of equations (21), (22) form a linear space. Now, define 
the scalar product 

(g,f)= fdfi(z,z)^)f(z); (25) 



then the entire functions f(z) of order p < 2 form the Fock space T [10,11], i.e. the Hilbert space 
of harmonic oscillator. We conclude: the phase functions f{z), describing non-equilibrium distri- 
butions of harmonic oscillator, can be identified with the probability amplitudes. For example, it 
is well known [10,11] that functions Z n (z/yh) give the orthonormal basis in the Hilbert space 

[Zni Z m ) = o nm , (26) 

and we demand 

oo 

(/,/) = E|Cn| 2 = l. (27) 

ra=0 

The latter equality means that |/(-2)| 2 may be identified with the probability density. Numbers 
\c n \ 2 also represent some (discrete) probability distribution. Evolution of a non-equilibrium dis- 
tribution (i.e. of \f(z) | 2 ) is given by evolution of a single function f(z). 

It is easy to see that smallness of c in (19) does not impose limitations on the set of functions 
e cz as a basis in T. It is known [10] that for arbitrary small e > vectors e CkZ , \ck\ < e form 
a (over) complete basis in JF if an infinite sequence of points c& in the complex plane converges 
to some limit Cq (one can omit from the sequence any finite number of points). It means that if 
some function ip(z) G T is orthogonal to all the functions e CkZ then ip(z) = 0. Indeed, functions 
e cz = (z\c) = c(z) play the role of (^-functions in J 7 , and (ck,ip) = (ck\ip) = i>(ck) = 0,VA;. But 
■ip(z) is entire function, so ip(z) =0. 

Commutation relations. It is well known that for the Fock space with measure (17) operator 
z defined as zg(z) = zg(z) has the Hermitian conjugate z = hd/dz (zf(z) = zf(z); it can be proved 
by integration by parts in (25) [10]). The commutation relation 

[z,z} = h (28) 

allows to identify z, z as the creation and annihilation operators for quantized harmonic oscillator. 
Then, the canonical variables q,p also become operators 

[q,p] = ^[z + z,z-z]=ih. (29) 

Remember, that we are working in the framework of pure classical theory. 

The Schroedinger equation. The classical equation of motion (21) can be identified with 
the Schroedinger equation if we multiply it by ih: 

ihZ n (z) = H C iZ n (z), H C i = ihH d = hud + d; (30) 

here a = d/dz, d + = z (because they act on functions f(z)). 

The spectrum of operator Hqi differs from that for quantized oscillator E n = hu(n+ 1/2). We 
would obtained correct answer using the Hamiltonian (20). Surprisingly, Eq. (20) automatically 
gives correct ordering of operators. It is important to elucidate the nature of the "quantum 
energy" hu/2. Functions Z n correspond to classical periodic motion with frequencies ion in the 
ID ip space (non-equilibrium states). But there is another canonical variable \z\, characterized 
by the Gibbs distribution, which also contributes to the total energy. So, the "zero energy" hoo/2 
should be attributed to the universal influence of the thermal bath. 
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Restoration of the equilibrium state. To model restoration of the Gibbs distribution one 
may introduce, for example, friction. The equation of motion for harmonic oscillator then reads 

q + aq + uj 2 q = 0, (31) 

where a > specifies friction. For infinitesimal a the general solution of (31) is 

q(t) = Cie -^-W2)t + C2e ^+W2)i ? (32) 

Evidently, there is only one equilibrium distribution. We see that the classical deterministic motion 
disappears in the limit t — ► oo (t > t r ~ a~ 1 ). It means that in this model of quantum oscillator 
all the probability amplitudes tend to zero when t — ► oo . 

5 Linear chain of harmonic oscillators in a heat bath 

It is not difficult to show that the linear chain of harmonic oscillators in a heat bath models 
standard ID quantum field theory. Again, it is the Gibbs distribution that is responsible for the 
non-trivial measure (see (17), (25)) in the Fock space. 

Hamiltonian of a chain of harmonic oscillators can be written in the form 

H=W(- + l(<ln-<ln-iY + iql) (7,7 >0). (33) 

z n m 

In normal coordinates a(k),a*(k) it reads 

1 ? 
H=- dku(k) {a*{k)a{k) + a(k)a*(k)) , (34) 

-A 



v^ y Jmu(k) 



where 



A A 1 

q n = dku(k)(p* n (k), p n = dkp{k)ip n {k), ip n {k) = — ==e^ k , A=-: (35) 

-A -A 

here a is the distance between the neighbour oscillators, u 2 (k) = 7/ra + 4(7/771) sin 2 (7r/c/2A). In 



7T 



the limit a — > 0, n — > 00, an — > x, a 2; y/m — > 1, 7/777 = M 2 , q n Jm/a — > (p(x,t), one obtains the 
ID theory of free scalar field tp with mass M. It is important that 

1) the Hamiltonian (34) presents a set of noninteracting oscillators with frequencies u(k), 

2) this theory is Lorentz invariant (the Lagrangian is (0 2 — ip' 2 — M 2 (p 2 )/2), so all the fields 
a(k),a*(k) with different k transform one into another by the Lorentz transformation. 

Evidently, the Gibbs measure ~ exp(— /3 J dkui(k)a*(k)a(k) differs from the Fock space measure 
exp(— / dka*(k)a(k)/h). But in fact the difference is formal. Introducing new variables a(k) = 
\0^.a(k), \ k = uj{k)/uj, uj = uj(0), one rewrites (34) 

A 

H = u I dkh*(k)a(k). (36) 

-A 



The Gibbs measure now takes the form 

TT \~ x 3ij — —i-Z e -1 ^ / dfc °*( fc )°( fc ) (37) 

k ' 2v ™ 

(the integration limits are omitted). The infinite factor T7 A/T 1 is n °t essential -- it is eliminated 
by normalization. The corresponding change of functionals $[0] — > $[A _1//2 a] leads to redefinition 
$[A~ 1//2 a] = $[a]. Thus, we come to the scalar product for complex valued functionals $i,$2 in 
the Fock space 

($ 1)$2 ) = /TT da *(*) * & ( fc ) e -/^ WaW /^ $ r a ] (3g) 

7 fc ^ 

(in the limits a -»• 0, A -»■ 00). From (38) it follows that a + {k)$[a*} = a*(k)$[a*], a(k)$[a*] = 
hS$[a*]/Sa*(k) : and [a(k), a+(k')] = h5(k - k'). 

This construction elucidates concepts of space, quantum mechanics and quantum fields. The 
physical space is nothing but the chain of harmonic oscillators, fields are excitations of the chain, 
particles (quanta) are one-particle excitations of the fields, and quantum theory appears as a 
method appropriate for description of non-equilibrium distributions characterizing the system. 
One should distinguish "the physical space" from the outer one in which it is embedded. The 
physical space is the linear chain of oscillators. Excitations of the oscillators ("matter") propagate 
in the physical space. The outer space is analogous to the multidimensional space (10D or 26D) 
in superstring theory. The 3D physical space also can be modeled in this way as a 3D network 
made of superstrings [2] . 

It is important that here the neighbours of an oscillator are responsible for taking it out from 
the equilibrium. The Planck constant H = l//3u is universal for all the oscillators, i.e. it is a 
fundamental constant of the theory. Distributions 

'^]^[a]e'I dka '' (k)a(k)/h (39) 

describe non-equilibrium states of the system (because they differ from the equilibrium one). 
To describe the evolution of density p[a*,a] = &[a]$[a] it is enough to describe the evolution 
of $[0], i.e. the evolution of non-equilibrium states is described by probability amplitudes (wave 
functionals) $[a]. In this model the Planck constant h and the wave functions (functionals) appear 
simultaneously. The chain of harmonic oscillators models also the bosonic string. 

6 Conclusion 

It was shown that small deviations from equilibrium states of harmonic oscillator in a thermal 
bath are described by probability amplitudes (in the time interval t r > t> uo^ 1 ). Thus, quantum 
mechanics can be modeled in framework of classical theory. Of course, it does not mean that one 
can obtain quantum theory of e.g. electron from classical theory of a pointlike particle. The result 
is meaningful only on the fundamental level of the Planck scales. For example, a chain of such 
oscillators in a thermal bath models one-dimensional quantum field theory (or quantized bosonic 
string). Excitations of the chain (quanta, particles) move along the chain — they model "matter". 
In this case, one simultaneously models both the ID physical space and quantum mechanics. 



Another important feature of the model — natural appearance of the Planck constant h (see 
(11), (17)). The probabilities, probability amplitudes and constant h appear as consequences of 
putting oscillator into a heat bath. The other attributes of QM (operators, complex functions) 
are indispensable features of CM. Actually, there is no need to introduce a thermal bath by 
hand. Complex ergodic systems are stochastic by themselves, i.e probabilities pervade Hamiltonian 
mechanics, and the Gibbs distribution also follows from CM. The "analogy of the evolution of 
classical correlation functions with quantum mechanics" was discussed in [12] (without introducing 
h; see also [13]). 

It follows from (32) that in the (1 + 1) universe, modeled by a chain of harmonic oscillators, 
quantum description will be actual only in the time interval (£ r ,u; _1 ); here t r is the relaxation 
time of the chain excitations (if it doesn't exceed that of the oscillators). When t — > oo the system 
tends to the equilibrium state ("vacuum"). Any coherent motion dies away. This is another kind 
of "death" of the "universe" in contrast with the "thermal death" of XlX-th century. In the 
latter case matter becomes thermalized, while here matter "disappears" because this is nothing 
but quantum excitation of the chain (i.e. vacuum) specified by some probability amplitudes. 
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